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Abstract 

We describe a unified approach to calculating the partition functions of a general 
multi-level system with a free Hamiltonian. Particularly, we present new results for 
parastatistical systems of any order in the second quantized approach. Anyonic- like 
systems are briefly discussed. 



Introduction. Statistics plays a fundamental role in determining macroscopic 
or thermodynamic properties of a quantum many-body system. So far there has 
been firm experimental evidence only for two statistics:Bose-Einstein and Fermi- 
Dirac.Much effort has been made to construct theories with generalized statistics 
and to search for deviations from Bose and Fermi statistics. Parastatistics [1-3] were 
the first consistent generalizations of Bose and Fermi statistics in any spacetime di- 
mension. They are invariant under the permutation group. The N-particle states can 
occur in a large class of representations (REP's) of the permutation group 5jv,i.e. 
not only in totally symmetric REP (for bosons) or totally antisymmetric REP (for 
fermions). Parastatistics can be related to theories obeying ordinary (Bose, Fermi) 
statistics and carrying suitable internal symmetries [3]. Hence a natural question 
arises whether there are systematic ways of distinguishing between the two classes 
of theories. 

Other generalized statistics have also been suggested, such as infinite quon statistics 
[4,5], a new version of parastatistics [6], fractional anyonic (braid) statistics [7,8] and 
other generalized statistics related to algebras of creation and annihilation operators 
with Fock-like representations [9-12]. 

To investigate physical consequences of such generalized systems, it is important to 
know the corresponding partition functions. There are few partial results for some 
parastatistical systems [13, 14]. Our motivation is to present a unified approach to 
calculating the partition function of a general multi-level system with a free Hamilto- 
nian. Particularly, we present new results for parastatistics and a correction to results 
for two-level systems for parafermions given in [13]. Anyonic-like operator algebras 
are briefly discussed. 
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The partition function for a general multi-level free system. The sta- 
tistical average of an observable, described by an operator O in a given ensemble, is 
defined as 

TrOe~^ H 1 , n aH 

<°>= -j-^r = z Tr(De ■ W 

Here Z is the thermodynamic partition function for a multi-level system described by 
M independent creation (annihilation) operators a\ (di),i = 1,2, ..M. The operator 
algebra is defined by a normally ordered expansion T (generally, no symmetries are 
assumed) [11]: 

^a] = r i:j (a\a), (2) 

with well-defined number operators iVj,i.e. [A^,a]] = aj^, [Ni,aj] = — and 
[Ni, Nj] = for i, j = 1, 2, ..M. We assume that there is a unique vacuum |0 > and 
the corresponding Fock-like representation. 

The scalar product (bilinear form) is uniquely defined by < 1 >= 1 , the vacuum 
condition Oj|0>=0, i = l,2,..,M and Eq.(2). A general N-particle state is a linear 
combination of monomial state vectors (a| x • • • a\ N \0 >), ii, • • • , ijv = 1, 2, ...M. 

We consider Fock representations with no state vector of negative squared norm 
but we might have norm zero vectors. These null-vectors will represent the only 
relations between the creation (annihilation) operators on the associated (non de- 
generate) quotient Fock space . 

For fixed sequence of indices i\---iN, we write its type as l ni 2 n2 ...M nM , where 
ni,n 2 , ■■■,riM are multiplicities of appearance 1,2, ...M in ii • ■ -i N . Note that rii > 
and YjiLi^i = iVThen, there are N\/ni\n 2 \...n M ^ (in principle) different states 
(aJi ' ' " a lj° >)labelled by all sequences of type l ni 2 n2 ...M nM and let A r (rn, ■ ■ -n M ) 
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be the Gram matrix of their scalar products. The number of linearly independent 
states among them is given by d^ ll .. = rank[A r (ni, ■ ■ -n M )], satisfying 

All the quantities nM completely characterize the statistics,the partition func- 
tion and the thermodynamic properties of the above free system realized on the 
corresponding Fock space. The free Hamiltonian is then 

M 

H = Y,E t N t , (4) 
i=i 

where is the energy of the i th level and N t are the number operators counting 
particles on the i th level. (Note that the statistics, i.e. the numbers d nij ... nM do not 
determine uniquely the whole operator algebra, Eq. (2).) The partition function for 
a free system described by the T algebra, Eq.(2),is given by 

oo 

Z T (x 1 ,...x M )=Tre-f ,Ho = £ Z T N (x u ...x M ) , (5) 

N=G 

where 

4(^1,-%)= E <,..n M ^---^ M - (6) 

niH hn M =N 

Here .. can be considered as the degeneracy of the state with the energy 
E = E,=i ,where ^ = e~^,(3 = ^. 

The central problem is how to compute d£ nM , ni, • • • , > O.The number of all 
allowed N-body states distributed over M energy levels is given by 

D r (M,N)= J2 <,...„ M ^(M^)- (7) 

nH hn M =JV M 

The partition function for permutation invariant multi-level systems. If 

the relations T in Eq.(2) are invariant under the permutation group SM,then the 



matrix A r (ni, ■ ■ ■ um) and its rank dj nM depend only on the collection of the 
multiplicities {ni, um}, which written in the descending order Ai > A 2 > ... > 
Am > , | A | = Y^iLi \ — N, give rise to a partition A of N,i.e. d^ M = <i A r and 
A T (n 1 , ■ ■ -n M ) = A{. 

If Ai = A 2 = • • • = Ajv = 1, Ajv + i = ••• = Am = 0, the correspponding Young 
diagram, denoted by 1^, is a column of N boxes. The corresponding N\ x iV! generic 
matrix is denoted by ^4^. All other matrices A^ , |A| = for any partition A of N 
are easily obtained from the generic matrix A\ N - By symmetry of T, it follows that 
A\ N can be written as 

A\s = £ cl R(n) , (8) 

where R is the right regular representation of the permutation group Sn and 
are complex numbers. If the operators a« and a\ are adjoint to each other, then 
the matrix A\ N is hermitian. Of course, d\ N = rank[A\ N ] < N\. The matrix 
A\ N commutes with every permutation matrix in the left regular representation. 
Hence the nondegenerate quotient Fock space splits into the sum of irreducible 
representations (IRREP's) of S N , and we can write 

^ = £n r (//)*>v, (9) 

where n r (/i) > is multiplicity and K^n (Kostka-Foulkes number) is the dimension 
of the IRREP fx of the group Sn- The sum in Eq.(9) runs over all partitions fx of 
N, i.e. all IRREP's. The multiplicities n r (/i) can be determined from the spectrum 
of the matrix A\ N - /,From Eq.(9) and Sm invariant structure of partition function 
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Zff(xi, ...xm) it follows that 

4 = E^)^,a, (io) 

with the same n r (/z) as in Eq.(9) and where are the Kostka-Foulkes numbers 
enumerating semistandard(column strict) tableaux of weight A and shape /z. Clearly, 
< K u,i N = dimfi.The set of all n r (/i) completely determine the statistics, 
the partition function and the thermodynamic properties of the free Sm invariant 
system. The Sm invariant partition function, Eq. (6), becomes 

Z^(x!,...x M ) = d{ J2 x i 1 '" x m 1= 2 d T x m x (x u ...x M ) 

A;|A|=JV 7r€5 M A;|A|=7V 

Z! n r (n)s li (x 1 ,...x M ), (11) 

H-M=N 

where m\(xi, ...%) is the monomial Sm symmetric function (summed over all dis- 
tinct permutations of (Ai, ...Am), and s /1 (xi, ...%) are the Schur's functions [15], sat- 
isfying 

s M (xi, ...x M ) = K i*,* m \( x i> - x m) (12) 

Note that for the numbers D T (M, N) of all allowable N-particle states of a permu- 
tation invariant M-level system we have: 



D r (M,N) = ]T d{m x (l,l,...l) = £ n r (/^(l,l,...l) 

A;|A|=7V ^ ^ ' W |/x|=iV ' ^ ' 

= Z n r (fi)dim{fi} (13) 
m;|mI=^ 

since s M (l, 1, ...1) = dim{fj,}, where {//} denotes the IRREP of S77 (M) corresponding 

M 

to the Young diagram {//}. 
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Remark. If we consider the interactions of particles described by the Hamiltonian 
H = H(Ni, ..N M ), where iV; are number operators,then 

Z F n(H) = E € u ... nM e-^ n —l (14) 

niH ^n M =N 

If the algebra r,Eq.(2),is invariant under the permutation group SV^then 

Z V N {H) = E rf A E e-W^W'-^wl (15) 

A;|A|=JV vr€S M 

Multi-level parastatistical systems. Here we analyze multi-level para-Bose 
and para- Fermi systems since they are important examples invariant under the per- 
mutation group S^-The operator algebra corresponding to parastatistics of order p 
[1-3] is defined by trilinear relations 

[ajcij ± aj-aj, aj] = (2/p)5 jk a\, Vi, j, k = l,2, M. (16) 

Its Fock representation satisfies the following conditions: 

< 0|0 >= 1 ai |0>=0 

a ia ]\0 >= 5ij\0 > i,j = 1, ...M (17) 

with |0 > denoting the vacuum state. The upper (lower) sign coresponds to para- 
Bose (para- Fermi) algebra, and p is the order of parastatistics. The Fock space does 
not contain any state with negative squared norm if p is a positive integer [3]. The 
consistency condition are the following [1,3]: 

[at,[at,4]±] = Vi,j,k, (18) 
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where the upper (lower) sign coresponds to parabosons (parafermions) . Note that 
Eq.(18) does not imply eq.(16). For p < N,in the N-particle space, there are ad- 
ditional null-states leading to relations not contained in Eq.(18).For the para-Bose 
and para-Fermi algebra,Eqs.(16),(17) can be presented in the form of Eq.(2), [11]. 
The matrices A p ^ e , \\\—N can be calculated recursively using the following relation 
[11,12]: 

N 

aid], ■■■ai N \o> = e <W fe_1 a h ■ ■ • "!* • • • a L 1° > 

k=l 

N fc-1 

- (2/p) E s m E e ' a l • • • 4 • • • 4-i4 a L+i • • • a L i° > • ( 19 ) 

where e = =Fl , the upper (lower) sign is for parabosons (parafermions), and the d 
denotes omission of the corresponding operator a. 

Case M = 2. Two-level parastatistical systems have already been studied in 
Ref.[13]. But the method in [13] is not completely correct (in fact the results 
for parafermions are false -see (20) bellow). For M = 2, from Eq.(16) it follows for 
parabosons (af) 2 * 3 ] = a \{ a \) 2 -,h3 — 1)2, and for parafermions (a|) 2 a] = 2a|a]a| — 
at(a}) 2 ,i,j = l,2. 

Then every N-particle state, N = ri\ +n,2,n 1 > n 2 , and ni,n 2 fixed, can be expressed 
in terms of the following n 2 + 1 vectors [13]: 

(air(4r |0 > ; (a\r- l aWMr- l \0 >;•••; {a\) n ^ a\a\ ■ ■ ■ a\a\ |0 > . 

S v ' 

2n 2 

In order to find the number c^' e (< A 2 + l)of linearly independent states of type A 
,we have studied the associate (A 2 + 1) x (A 2 + 1) Gram matrices up to iV = 6. For 
example, for iV = 6 there are three cases, i.e. A 2 = 1,2,3 with 2x2,3x3,4x4 types 
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of matrices, respectively. We found that for parabosons (with p > 2),d\ — A2 + 1, i.e. 
the result does not depend on p and coincides with the number of allowed IRREP's 
\x for fixed Ai, A2 in the decomposition, Eq.(lO). However, for parafermions p > 1 we 
found up to N < 6 (in contrast to [13]): 

f ifp< Ai 

= I p - Ai + 1 if A x < p < Ai + A 2 = N (20) 
( A 2 + l if p > Ai + A 2 = N 

i.e. the result does depend on p if p < N,in contrast to [13] and just coincides with 

the number of allowed IRREP's ji for fixed Ai, A 2 and Ai < p in the decomposition 

Eq.(10). Furthermore, if number of parts in A (the lenght of A), /(A) < 2, then 

K[i\ = 1 for any of the 2A 2 + 1 allowed /i, i.e. for /i = (Ai + A 2 , 0), •••,// = (Ai, A 2 ), 

[15]. This means that every IRREP ji for A, with /(A) < 2, is filled only with one 

state. Hence, d p ^ e = X^^ p ' e (/-0 and sum runs over allowed IRREP's for parabosons 

or parafermions. 

The results for multiplicities n p ' e (fi) of IRREP fi in d p ^ € Eq.(10), (checked on a 
computer up to iV < 6, /(A) < 2) are for parabosons p > 1: 

n vB, n = f 1 if < min(2,p) 
I otherwise 

and for parafermions p > 1: 

n pF/ n = f 1 if < 2 a ™ rf ^(/^ T ) < P (22) 
I otherwise 

where /z T is the transposed tableau of /i. Note that the results, Eqs. (20-22), for a 
given A = (Ai, A 2 ) are valid for all M > 2. 

Starting with Green's second quantized approach eqs.(16)-(18) and using the theo- 
rem proved in next subsection, we obtain the following : 



General result :n p,e (/i) = 1 for all allowed IRREP's /i in the decomposition of c^' e ,eq.(10). 

In particular for p > N = = |A|, the multiplicities and dimensions for parabosons 

and parafermions coincide, i.e. n pB (n) = n pF (fi) and d p x B = d p x '. 

For the partition function for two-level (M = 2) parabosons of order p > 2 we then 

have: 

Z pB (x u x 2 ;p) = 1 -. (23) 

(1 - Xi)(l - x 2 )[l - XiX 2 ) 

i.e. it is independent of p (in agreement with [13]). The partition function for two- 
level (M = 2) parafermions of order p satisfies the recurrence relation, which follows 
from ^' + , Eq.(20): 

ZP F (x 1 ,x 2 ;p) = Z pF (x 1 ;p)ZP F (x 2 ;p)+x 1 -x 2 Z pF (x 1 ,x 2 ;p-2), (24) 

where 

p+i _ 1 

Z pF [x-p) = p>0. 

x — 1 

The solution of eq. (24) is given by 

[f] 

Z pF ( Xl ,x 2 ;p) = Y.^2) k Z pF {x l -p-2k)Z pF {x 2] p-2k) 

k=0 

= 1 r l-(W +2 X\ +2 ~XI +2 

(1 - Xi)(l — X 2 ) 1 - XiX 2 Xi~X 2 

This result is new and differs from that in [13]. In particular, we find for p — 1,2: 

Z pF (x 1 ,x 2 ;l) = (l+x 1 )(l + x 2 ), 

Z pF (x u x 2 ; 2) = (1 + x ± + x\){\ + x 2 + x\) + x ± x 2 , 

and for p > 3, Z pF (xi,x 2 ;p) can be easily found using Eq.(24) or Eq.(25). In the 
limit E 2 — > oo, x 2 — > 0,Eq.(23),(25) reduces to the partition function for one-level 
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parabosons (parafermions) . Up >> 1, then Z pF (xi, x 2 ; p) ~ Z pB (xi, x 2 ; p\ However, 
for p = oo the para-Bose (para-Fermi) algebra (16) becomes Fermi (Bose) and hence 



Generally, there is no smooth transition from generalized statistics to Bose or Fermi 
statistics, or from "higher" to "lower" statistics. 

General case M > 2. Here we generalize the results of Ref.[13] to any multi-level 
system and the results of Ref . [14] to parastatistics of any order following the Green's 
second quantized approach, eqs.(16)-(18). 

Theorem. Let us consider the iV-particle states of type l n ^2 n2 ■ ■ ■ M nM , £*=i n { = N. 
If the order of Green's parastatistics defined by eqs.(16)-(18) is p > N, then the 
number d\ of linearly independent physical states is 



where denotes the Kostka numbers and A is partition of N with parts rii, n 2 , n M 
Sketch of proof. Let us consider all possible states of the type l ni 2 n2 ■ ■ ■ M" M , 
HiLi n i — N. Our aim is to write the canonical basis for this subspace which gen- 
eralize the corresponding basis for M — 2, Ref. [13]. 

We shall transform any given TV-state by using the para-Bose (para-Fermi) rela- 
tions, eq.(18) in the following form (as reduction rules in which, for simplicity, we 
abbreviate a\ to k etc.): 




Z pB (xi,x 2 ;p = oo) = (1 +xi)(l + x 2 ) = Z F (x 1 ,x 2 ). 



(26) 



kji 



ikj + e jki — e ijk 



kki 



e ikk + (1 + e) kik, 1 < i, j < k 



(27) 
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(here e = +/— refers to parafermions/parabosons). We first start shifting the 
rightmost index M to the right. This will end with M either in the last or next to 
last right position * M or * Mi, i < M. Then we apply the same procedure to the 
second rightmost index M by shifting it to the right as before. After moving all of 
M's, we end up with states like * M" M , * M i x M nM_1 , • • •, * M i x M i 2 ■ ■ ■ M i nu , 
i.e. with riM + 1 types of states. 

By applying the approach of Bergman, Ref.[16] to the relations (27), we find new 
important relations (no further relations except these exist if p > N) 

kjki = kikj + jkik — ikjk, k > j > i, (28) 

which together with (27) ensure the diamond condition for the para-Fermi/para- 
Bose algebras. By applying (28) to the above types of states we further reduce them 
to the states satisfying 1 < i\ < i 2 < • • • < i nM < M. 

After exhausting M's we perform the same procedure on (M — l)'s , (M — 2)'s 
,.. and finally on 2's. In this way we obtain a canonical basis which generalizes 
the basis for two-level system of the previous subsection. We point out that the 
dimensions d\ satisfy the same recursion relation as coefficients a\, ( a x = ^« K^x) 
in the expresion fl^i 1/(1 — x %) \~\i<j 1/(1 — = a \ m x- Hence, we conclude 
that the eq.(26) in the theorem is proved. (Alternatively, one obtains one to one 
correspondence between the states in canonical basis and terms in the expansion of 
E M s IM (x 1 ,...,x M )-) 

Comparing the eq.(26) with eq.(10), we obtain n{ji) = 1 for all allowed IRREP's /i, 
i.e. /(/i) < M for para-Bose and para- Fermi case. 

Furthermore, if p < N then the Green's parastatistics imply that only IRREP's 
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with < p for para-Bose and /(/i T ) < p for para-Fermi are allowed. These are 
additional restrictions on the multiplicities n(/i) in eq.(10). We point out that these 
restrictions are not put by hand, but they are intrinsicaly contained in eq.(16). The 
complete proof of the theorem and corrolaries will be presented elsewhere. 
This theorem enables a simple generalization of Eq.(23).The partition function for 
multi-level parabosons for p > M is: 

M i M 1 

Z? B ( Xl , ...xu-p >M) = J2 s x (xi, -..%) = II 77— TT II n _ TrV ( 29 ) 

A i=l V 1 X V i<j \ L x t x j) 

The identity used in Eq.(29) can be found in [15]. The results for parastatistical 
systems of order p = 2 [14] can be easily obtained from Eqs. (10,11) and our eq.(26). 
We point out that the parastatistics in the first quantized approach of Ref.[17], 
[18] leads to the same partition functions as we have found. However, these two 
parastatistics basicaly differ: (i) since within the first quantized approach [17], [18], 
the Green's trilinear commutation relations , eq.(18) cannot be obtained, and (ii) 
the probabilites of finding IRREP T c in r a <g> T b = Y, c T c and finding T a <g> T b in 
T c = J2( a ,b) r a <E> r b calculated within the approach of [17], [18] would be completely 
different from those calculated in the second quantized approach. Hence these two 
approaches are not equivalent. The complete analysis will be published separately. 
Anyonic — likemulti — level systems. Here we study another class of multi-level 
systems described by algebras with Fock representations,which possess anyonic-like 
statistics and are not invariant under the permutation group. We do not assume 
any symmetry requirements . Instead, we assume for % ^ j: 

aid] = e^'ajai, faj = -</>ji e R, j, i,j = 1, 2...M (30) 

and aj, a\ are adjoint to each other. It is easy to show that d nit ... nM = 1 for any set 
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of mutually different indices i±, ■ ■ - In- In this case, any monomial state consisting of 
permuted indices - • - in is ,up to unit phase, equal to (a^ • • • a\ N |0 >). Therefore, 
we find the anyonic-like commutation relations [5]: 

a\a] = e~ l(pij a]al, = e~ l(p ' j ajcii i^j- (31) 

We have assumed no specific relation a^af = T(a\a) for the same indices. The 
independent number operators JVj are assumed. Then for a multi-level anyonic 
system with the free Hamiltonian H = Y^Li EiNi, the partition function is 

M 

Z(x u ...x M ] <t>ij) = II z ( x *>Pi)> ( 32 ) 
i=i 

where Z(x;p) = (1 — x p+1 )/(l — x) is the one-level partition function of a single 
(parafermion-like) oscillator satisfying a p+1 = 0, p G N. Such a system is physically 
equivalent to the set of M commuting ( or anticommuting ) single-mode generalized 
oscillators a\, fij. The two sets of operators {c^} and are connected by a 
generalized Jordan- Wigner transformation [19]. 

Starting from Haldane's definition of generalized exclusion statistics [20], Karabali 
and Nair [21], using additional assumptions derived an operator algebra of the type 
given by eq.(30) with <pij = -^isign{i — j), i ^ j and af +1 = 0, for Vi From our 
results it follows that the partition function corresponding to the system described 
by the Karabali-Nair algebra with the free Hamiltonian H = J2fLi Ei iV; is 

M 1 _ t p +1 

Z(x ll ...x M ;p) = H-—^—. (33) 
f=i 1 ~ x i 

The system with this partition function was analyzed in [22]. However it does not 
reproduce the Haldane-Wu interpolating distribution. 

13 



Finally,let us mention algebras with no independent number operators N i} but with a 
well-defined total number operator N , [N, a\] = a], [N, a*] = — a« for 2, j = 1, 2, ..M. 
Let us consider a free Hamiltonian of the form H = EN. Then the partition 
function is 

00 

Z r ( Xl , ...x M ) = £ ^ r (M, JV)e-^, (34) 

Af=0 

where D r (M,N) is number of independent states (a\ x ■ ■ -a\ N \Q >), 
*i, • - - , «jv = 1,...,M. If a^- = aj-Oi, then < D(M,N) < D B (M, N). Such 
an example is the system derived from the Calogero-Sutherland model analyzed in 
[23]. A complete analysis including thermodynamic properties of this and other 
generalized multi-level systems will be treated separately 

Acknowledgements. The authors thank I. Dadic and M. Milekovic for useful 
discussions. 
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